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^, ABSTRACT 

o : 
o . 

CN ■ Fully relativistic and causal equations for the flow of charge in curved spacetime 

^ ■ are derived. It is believed that this is the first set of equations to be published that cor- 

Q ■ rectly describes the flow of charge, and evolution of the electromagnetic field, in highly 

^ ■ dynamical relativistic environments on time scales much shorter than the collapse time 

{GM/c'). The equations, therefore, will be important for correctly investigating prob- 
>- ■ lems such as the dynamical collapse of magnetized stellar cores to black holes and the 

^ ■ production of jets. Both are potentially important problems in the study of gamma-ray 

O ' burst engine models and in predicting the dynamical morphology of the collapse and 

^ ■ the character of the gravitational waves generated. 

g ■ This system of equations, given the name of 'charge dynamics', is analogous to 

■ those of hydrodynamics (which describe the flow of mass in spacetime rather than the 

p flow of charge). The most important equation in the system is the relativistic general- 

Q ■ ized Ohm's law, which is used to compute time-dependent four-current. Unlike other 

+i . equations for the current that are presently in use, this one ensures that charge drift 

. velocities remain less than the speed of light, takes into account the finite current rise 

> . time, is expressed in a covariant form, and is suitable for general relativistic compu- 

^ . tations in an arbitrary metric. It includes the standard known effects (Lorentz force, 

. Hall effect, pressure effect, and resistivity) and reduces to known forms of Ohm's law 

in the appropriate limits. In addition, the plasma particles are allowed to have highly 
relativistic drift velocities, resulting in an implicit equation for the 'current beaming 
factor' Yg. It is proposed that, short of solving the multi-fluid plasma equations or the 
relativistic Boltzmann equation itself, these are the most general expressions for rela- 
tivistic current flow in the one-fluid approximation, and they should be made part of 
the general set of equations that are solved in extreme black hole accretion and fully 
general numerical relativistic collapse simulations. 



Subject headings: magnetohydrodynamics — black holes — relativity — gravitational 
collapse — gamma-ray bursts 
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1. Introduction 



This paper is another in a series whose goal is to establish the mathematical, physical, and nu- 
merical tools necessary to understand and simulate the formation of black holes and the production 
(through electrodynamic processes) of relativistic jets during that collapse. The full understand- 
ing of how black holes are formed, and how jets and gravitational waves might be generated in 
accretion and gravitational collapse, is currently one of the most challenging, and far-reaching, as- 
trophysical problems. It has important observational consequences for gravitational wave sources, 
gamma-ray bursts (GRBs), quasars and microquasars. Its solution will involve nearly every branch 
of theoretical astrophysics (nuclear and particle physics, electromagnetics, gravity [numerical rel- 
ativity], plasma flow, radiation transport and dynamics). And the formulation of the problem will 
require that these processes be expressed in a general relativistic framework that respects the princi- 
ples of causality and covariance on time scales considerably shorter than the light crossing time of 
the forming black hole (< < GM /c^). While there is a good understanding of how charge behaves 
near black holes in equilibrium situations, i.e. on times >> GM/c^ (Wald 1974; Lee et al. 2001), 
to date there appears to be no thorough treatment of charge flow and field evolution in strong grav- 
ity on very short time scales. Such a treatment should take into account the fact that the current 
rise time can be long compared to the characteristic time scale, properly account for particle veloc- 
ities (bulk, drift, and thermal) that can approach c, be expressible in an arbitrary spacetime metric, 
and also reproduce, in the appropriate limits, the standard known effects (including the Hall and 
pressure effects, not just electric acceleration, Faraday induction and Ohmic resistivity). 

The purpose of this paper is to return to the basic equations of general relativistic statistical 
mechanics and properly derive the covariant equations for the evolution of the current. Section 2 
discusses the framework of the problem and the current lack of a good equation for the current in 
highly relativistic situations. Section 3 sets up the relativistic Boltzmann problem and section 4 de- 
rives the relativistic plasma equations. Section 5 discusses relativistic one-fluid plasma theory and 
derives the generalized Ohm's law and charge dynamics. Section 6 shows that the theory reduces 
to various previous generalized Ohm's laws in the appropriate limits and discusses the applicability 
of the one-fluid theory. The most important results of this paper are the basic equations of charge 
dynamics (56 - 58) and (63) or, in component form, equations (68 - 71). 



2. Causal and Covariance Problems in Present Treatments of the Current 

The two governing sets of equations for the general numerical relativity problem of electro- 
magnetic black hole formation are the Einstein equations 

G = (1) 
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for the gravitational field and the Maxwell equations 

4% 

V-F = —J (2) 

c 

V-M = (3) 

(Faraday's and Ampere's laws) for the electromagnetic field. V = = d/dx^ is the four-gradient 
operator; G is the second-rank, symmetric Einstein tensor that describes the second derivatives 
(curvature) of the metric g; and F is the second-rank, anti-symmetnc Faraday tensor that describes 
the electromagnetic field. T is the energy-momentum-stress tensor that serves as the source for 

gravity and J is the four-current vector that serves as the source of the electromagnetic field. The 
Maxwell tensor M =*F (the dual of the Faraday tensor), has the same information as the Faraday 
tensor and also is anti-symmetric, but has the roles of electric and magnetic field reversed. Because 
of their geometric properties, these tensors satisfy the following Bianchi identities 

V-G = (4) 
V-(V-F) = (5) 
V-(V-M) = (6) 

which, from equations (1) - (3), give rise to the conservation laws of energy and momentum 

V-T = (7) 

and of charge 

V-J = (8) 



The solution of equations (7) gives the distribution of temperature T and four-velocity U(x), 
which is constrained to always have the absolute value of the speed of light 

U-U = -c^ 

The addition of particle/rest-mass conservation 

V-p^U = (9) 

where Pm is the rest mass density in the fluid frame, allows a solution for pm(x) to be found as well. 
Together r(x) and Pm(x) can be used to compute the state variables (pressure, internal energy, etc.) 
that close the energy and momentum conservation laws. 

The solution to equation (8), however, which can be written as 



V-(p,U+j) = 
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gives only one quantity, the charge density in the fluid frame, in terms of the spatial charge 
current 



P-J (10) 



where 



P = lu®U + g (11) 

is the spatial projection tensor orthogonal to the four-velocity unit vector eu = U/c. The symbol 
(g) signifies the outer tensor product. While j is a four- vector, it is constrained to have only three 
independent components by its orthogonality to U 

U-j = (12) 

However, none of the above equations can be used to specify these three components of the current. 

Previous treatments of black hole electrodynamics have made a variety of assumptions to 
determine j and thereby close the electromagnetic equations. One popular technique is to assume 
that the electromagnetic field dominates the dynamics and is time-independent, leading to the 
"force-free" condition 

J-F = (13) 

This assumption has been criticized (Punsly 2003) as leading to effects that violate causality. That 
is, the force-free condition is acausal and, therefore, not relativistically acceptable for forming 
black holes and GRBs. 

Another approach has been to assume general relativistic magnetohydrodynamics, as outlined 
by many in the past (Lichnerowicz 1967; Eckart 1940; Anandan 1984; Blackman & Field 1993). 
This approach relates the current to the electromagnetic field through a simple form of Ohm's law 

Tij = --F (14) 

c 

where ri(x) is the resistivity distribution of the plasma. In the nonrelativistic limit, this takes on the 
familiar form 

o( £ + - xB 



where a = 1 /r| is the plasma conductivity and /, E, V , and B are the current, electric field, velocity, 
and magnetic field three-vectors. 
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Ideal relativistic MHD is a further simplification current flow that is useful for highly con- 
ducting plasmas, such as those in most astrophysical situations. With r| — > 0, Ohm's law reduces 
to 

U-F = (15) 
which, in the nonrelativistic limit, becomes 

E = --xB 
c 

The ideal MHD condition is not so much an equation for the current as a condition on components 
of the Faraday tensor (the electric field): when the conductivity is high, the local electric field in 
the plasma shorts out, leaving only the EMF due to charged plasma motion in the magnetic field. 
The current itself is determined by first determining the electromagnetic field from equation (3) 
and then inverting equation (2) for J. Equations (13) and (15) appear very similar. They both 
state that the Faraday tensor is orthogonal to a four- velocity, either the particle drift velocity or the 
average particle velocity. However, they generate very different physics. 

Unfortunately, a criticism that is equally as harsh as that of the force-free condition can be 
leveled against the MHD condition, even the resistive version of it (14). The latter states that the 
application of an electromagnetic field instantaneously generates a current. However, this also 
is acausal. There is no immediacy in relativistic dynamics. A current builds up after a finite, 
albeit short, rise time. And, in relativistic flow, the time-dependent increase in a current might be 
interrupted by any number of other rapid phenomena, resulting in possible charge separation and 
time-dependent charge dynamics. 

A nonrelativistic version of the "generalized Ohm's law" is often used in laboratory plasma 
physics and is given by the expression (Rossi & Olbert 1970; Krall & Trivelpiece 1973) 

+ W-(yj+JV-PqVV) + Wpq 

= e(E+^xB + h'^xB - r[J^ (16) 
where V is the three-space gradient operator and the charge-weighted pressure (per unit mass) is 

Pq = Y^laPa/ma 

a 

qa is that particle's charge, is its mass, and pa is that species' partial pressure. The Lorentz, 
Hall, and resistivity coefficients are 

a 
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^ a 

Ti = v/e 



'-a 



where Ua is the number density of particle species 'a'. The resistivity term results from integrating 
particle collisions over velocity and approximating the result as an effective collision frequency v 



Equation (16) has most of the effects we are looking for in a description of charge dynamics (finite 
current rise time, Lorentz force, Hall effect, pressure effect, and resistivity). However, it is valid 
only to linear order inV/c (still acausal), only in the laboratory coordinate system (not covariant), 
and only in flat space. 

Ardavan (Ardavan 1976) derived a relativistic form of Ohm's law for a cold plasma (vanishing 
Pa and pq) in flat space. This expression turns out to have some errors, but is still useful for check- 
ing our results in section 5. (Expressions for a relativistic cold pair plasma also have been derived 
(Gedalin 1996; Melatos & Melrose 1996), but the Ardavan expression is a little more general and 
therefore more useful to us.) In our notation the corrected Ardavan equation for the current is 



where U and J are the temporal components of the four-velocity and four-current, respectively; 

and are the spatial three-vector components of those four-vectors (i.e., = yV = U^V/c). 
The errors that have been corrected are all on the right-hand side: a sign error in the Lorentz term 
and the addition of —hpqU^E /c in the Hall term. In the limit of nonrelativistic flow, — > c, 
— > PqC, V, and — > /. So, with these corrections, equation (17) reduces to (16) when the 
charge- weighted pressure pq = 0. Equation (17) is both causal and covariant, but it is valid only 
for Lorentz systems in flat space and only when the plasma is truly cold. It does not include effects 
that occur when the plasma has a relativistic temperature or a relativistic current drift velocity. 

The goal of this paper is to derive a description of charge flow that is valid in all relativistic 
situations — relativistic bulk flow, hot plasma, relativistic current drift velocities — and is causal, 
covariant, and valid in any spacetime metric. 




^[U0 3j + (J' 



e!.^[U^E + ^VxB] + [{f-pq\J^)E + ^JxB] - Ti^jj (17) 
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3. General Reladvistic Statistical Mechanics 

3.1. Phase Space and Particle Density 

Phase space Q. = Q.x <^ Q.u inherently eight-dimensional, not six. The generalized coor- 
dinates are the position and four- velocity (x^^jU^), where //,v = 0, 1,2,3. Q.x is a general curved 
spacetime with a global time coordinate = ct and three spatial coordinates. The volume element 

(where is the determinant of the spacetime metric) is an invariant over the entire spacetime. 
Therefore, in a general metric, dQ.x cannot be separated into globally invariant temporal and spatial 
parts. However, in each local Lorentz frame at a given point in spacetime, it can be separated as 



d^x = dx dTx 



where 



dx = JxV(u7c) 
dXx = {\x^/c)dx^dx^dx^ 

Each of the factors in equation (18) are Lorentz invariant (Misner et al. 1973). We therefore can 
define a Lorentz invariant three-space density of particles of species 'a' and four- velocity u in each 
local Lorentz frame 

^ (18) 
dYxdn.u 

in the ^even-dimensional phase space Tx®^u- 

While Clu is a four-dimensional velocity space, all particles are constrained to move on a 
three-dimensional hypersurface within that space (the "mass hyperboloid") T„, defined by the 
normalization of particle velocity 

u-u = -c^ (19) 

In regions of Q.^ outside of T[, the particle density vanishes, as there are no particles with u - u ^ 
—(?■ . Some treatments of relativistic statistical mechanics incorporate the constraint (19) into the 
Boltzmann equation directly. Indeed, one can separate d^u into temporal and spatial parts, in a 
manner similar to dtlx 

dilu = dzdTu (20) 



-8- 



where 

dz = (u7c)Ju" = J(u")V2c (21) 
dJu = dv}d\x^d\y' /{xx^/c) (22) 

Not only are de and dT^ Lorentz invariant, the product 

dJx dXu = dx^ d\^ dx^ du^ du^ du-^ 

is also Lorentz invariant (Misner et al. 1973). One therefore could define an invariant density 
in 5/;c-dimensional phase space instead of it^. However, in this paper it has been found to be 
more useful to use the seven-dimensional phase space and perform the velocity integrals over Clu. 
The mass hyperboloid then is enforced only at the end of the computation when the integrals are 
evaluated. This is accomplished by using a delta function to describe the lack of particles outside 
of the mass shell in the distribution K^. (See Appendix A). 

A final property of the particle density to note is that, as |u'| ^ oo (where /= 1,2,3), 
approaches faster than any power of u'. Therefore, )^a, multiplied by any power of u', vanishes 
on the hypersurface aQ.u (the boundary of Q.u)- 



3.2. The Relativistic Boltzmann Equation 

The density for each particle species obeys the relativistic Boltzmann equation 

d)^a 



dx 



= ^a,col (23) 



where T again is the proper time for particles in that region of phase space, and K a,coi is the number 
density of collisions per unit time of that particle species at that point in phase space (summed over 
all other particles of all other species) 



^a,col = -Y, [ ,{a-a')-^u^'ab{^,n,x',u')da' 



(24) 



with a being the particle acceleration caused by body forces (i.e., forces other than particle colli- 
sions). Equation (23) can be rewritten in generalized coordinates as 

S ^^j. — t^a,col 

where s'' = (x^, u^) with r = 0, 1, 6, 7 or, in geometric form, 

a.col 

(25) 
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and the gradient operators are 




For the electromagnetic field a is the Lorentz acceleration 

a = ^u-F 

Some treatments of general relativistic statistical mechanics (Andreasson 2002) also include the 
gravitational "force" by including Christoffel symbols in the acceleration. However, this is not 
necessary, nor really desirable, as we can implicitly take these effects into account by using only the 
geometrical form of the gradient operator throughout the derivation and then using the equivalence 
principal at the end to convert the equations to component form. As a body force, gravity is 
automatically included in the structure of whatever spacetime in which S/x is evaluated. 



4. The Multi-Fluid Equations 

The equations of plasma dynamics are generated by taking velocity moments of the relativis- 
tic Boltzmann equation, which gives rise to many hydrodynamic and thermodynamic quantities. 
Taking the moment involves multiplying equation (25) by a power of the velocity coordinate vector 
u and integrating over all velocity space Ci^. 



4.1. The Zeroth Moment: Conservation of Particles 

Multiplying equation (25) by unity and integrating over Clu produces the zeroth moment. One 
can show that the integral of the second term in that equation vanishes by first integrating by parts 



Jq.u ./f2„ Jsiu 



(26) 



Gauss' law then can be used to show that the first term in (26) vanishes because vanishes on 
the boundary 
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where <iE„ is the three-volume element on dClu. The kernel V„ • a in the second term in (26) also 
vanishes because the Faraday tensor is antisymmetric and independent of the velocity coordinate^ 

v.-i'-^) = 

Similarly the velocity integral of the right-hand side of equation (25) vanishes because ^a.coi is 
in a similar form to the a ■ V„ term in equation (25). (See equation 24.) Because x and u are 
independent generalized coordinates, V;^ • u = 0, so the velocity integral of the zeroth moment of 
equation (25) becomes simply 



idQu = (27) 
In order to interpret this equation, we need a coordinate gauge in which to express u. 



4.2. Velocity Decomposition and the Velocity Coordinate Gauge 

We choose to decompose the velocity coordinate u into four components, one along the 
center-of-rest-mass average particle velocity U 

U ^ La'naJan)!^addu ^^^^ 

(the justification for this choice of average velocity will be given later) and one orthogonal to U 
(the drift four-velocity v), yielding 

u = Y(U + v) (29) 

where the Lorentz factor is 

Y = -4(U-u) (30) 
and the relative spatial velocity coordinate is 

V = (P-u)/Y (31) 



^In component notation, 
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Note that v is still a four-vector, but it is constrained to have only three independent components 
by its orthogonality with U 

U-v = (32) 

It therefore is a velocity coordinate that spans Y„. In the rest frame of the fluid, the velocity 
components are 

u = (yc yv\ Yv^ yv^) (33) 
V = (0, w\ v^ v^) (34) 

and the constraint (equation 19) on the particle velocity u becomes 
as expected. 

With the velocity decomposition in equation (29), the equation of continuity (27) for particle 
species 'a' becomes 



where 



is the particle density and 



V^-n«(U + V„) = (35) 



fla = f y^adClu (36) 



ya = - i yy^adCiu (37) 



is the average particle drift velocity for species 'a' . The equation of continuity (35) is important for 
the conservation laws of rest mass and charge in one-fluid dynamics. Despite the Lorentz factor in 
the above velocity integrals, these quantities are, in fact, the familiar three-momentum integrals of 
standard thermodynamics, in which the single factor of y does not appear. (See Appendix A). 



4.3. The First Moment: Conservation of Particle Energy-Momentum 

The first moment of equation (25) generates a vector equation 

V;,-(i?„u®u) +u|^^-V„k4 = U^a,col (38) 
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which also can be integrated over dClu to yield 

V;,- [n'«U®U + naU®V'a + naV'a®U + nJ = — Ja-F - vna(U + Va) (39) 
'- niaC 

where 

n'a = f y'^adCl^ (40) 

v'„ = - / fy^adCiu (41) 

In deriving equation (39) we have substituted equation (29) for u in the first term of (38) and 
discarded a vanishing boundary integral that results from integrating the second term by parts. 

Note the extra factor of y in the integrals in equations (40) and (41) compared with (36) and 
(37). These are beamed quantities that give rise to relativistic internal energy, pressure, etc. The 
partial current is 

Ja = qa f U^ad^ = qana{U + \a) (42) 
which has components qafia (c, V^, V^, V^) in the fluid rest frame. A four- vector 

ja = qana^a (43) 

can be used to describe the U-orthogonal part of the current (cf., equation 10), with U • = 0- The 
partial pressure tensor (per unit mass) of species 'a' is 



Ua = f (P-v)®(P-v)i^ A 



= i i'{y®y)^adCiu (44) 

Note also in equation (39) that the collision term has been simplified to a collision frequency v 
times the integrated particle flux. 



5. The One-Fluid Equations 

At present, astrophysical simulation codes deal almost exclusively with the one-fluid equa- 
tions when dynamics are concerned. That is, the individual particle species dynamical equations 
are summed and solved as a single set of equations. Of course, many stellar evolution and collapse 
codes track the composition for different species of particle, but this is done usually for equation of 
state and composition purposes only, not to determine relative drift velocities of different charged 
species, for example. Therefore, for at least the near future, it still will be important to construct 
one-fluid equations for studies of processes like black hole formation and accretion. 



5.1. The Example of Hydrodynamics 



Before deriving the equations of relativistic charge dynamics, it is important to review the 
derivation of the equations of hydrodynamics in the presence of an electromagnetic field. This is 
more than just an exercise in "re-inventing the wheel". Reviewing this derivation will allow us to 
check that our equations and procedures are correct and define some quantities that will be needed 
later in the discussion of charge dynamics, and it will assist us in understanding the new set of 
equations in terms of the familiar hydrodynamic ones. The procedure is to simply weight the five 
equations (35) and (39) by the particle rest mass and sum over all species 

V,-p„U = (45) 

Vx-|[pm + e/c^]U(g)U + ^[U(8)H + H(8)U] + ppj = ij-F (46) 

where 

Pm = Y^nama (47) 

a 

is the total rest mass density, 

8 = Y^^a = Y,{na-na)maC^ = Y^maC^ i y{y-l)>^adClu (48) 

a a a •'^u 

is the internal (kinetic) energy density, 

H = J^namaC^V'a (49) 

a 

is the heat flux (including relativistic enhancement), and 

P ^ l.Pa^l,'^ iy'(yy)^adCiu (50) 

is the scalar pressure for an isotropic distribution in T„. Equation (45) is correct only if the rest- 
mass-centered drift velocity is zero 

J^/lamaVa = 

a 

and that is the case only if the average velocity is defined according to equation (28), thereby 
justifying our rest-mass-centered choice for U. 



-14- 



We can show that equation (46) is equivalent to equation (7) if we define the following energy- 
momentum-stress tensor^ 



T = T^^ + T^^ (51) 



1 



= [p^ + (e + /7)/c^]U®U + ^[U^H + H®U] + pg (52) 



rjiEM 



1 

47t 



F-F- l(F:F)g 



(53) 



and recognize from equation (2) that 



ij.F = -V^-T^^ (54) 

c 

The multi-fluid equations of section 4, therefore, reproduce the familiar equations of general rela- 
tivistic hydrodynamics. 

Note that the energy-momentum-stress tensor in equation (52) is the one for an ideal gas with 
heat flow. It does not contain terms for viscous momentum and energy transport, however. These 
latter terms would arise if we performed a more sophisticated treatment of the collision term. 
Also, unlike nonrelativistic treatments (which require the second moment for the energy equation), 
the treatment here derives the conservation of energy using only the first moment. This occurs 
because the equations are relativistic and use the four- vector u in the first moment instead of the 
three-velocity v. The conservation of energy equation can be extracted by taking the component 
of equation (46) along the average velocity^ 

l]-V,e + {e + p)V,-lJ = -j-F-U- V;,-H-H-(U-V;,U)/c2 (55) 

which is also known as the "first law of thermodynamics": the change in internal energy £ plus 
mechanical work is given by Ohmic heating minus losses due to heat conduction. 



^In component notation the inner and scalar products of two tensors are 

[F.G]«P = F«^gPv 
F:G = F^G^ 



^The Euler equations also can be extracted by projecting equation (46) with the projection tensor P, but they are of 
no interest in this paper. 
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5.2. Charge Dynamics in Geometric form 

A similar set of equations to (45) and (46) can be generated by weighting (35) and (39) with 
the particle charge qa rather than the rest mass. The results are the equations of charge dynamics: 
the conservation of charge 

V;c-J = V^-(p^U+j) = (56) 
(equivalent to equation 8) and the relativistic generalized Ohm's law 

1 



(57) 



■(U + /ij)-F-Ti(p,U+j) 

where the charge-current-pressure tensor is given by 

C = [p^ + (8^ + p^)/c2]U0U + U®j'+j'®U + p5g (58) 

Note the appearance of a beamed current}' in the charge-current-pressure tensor while the source 
terms involve the unbeamed current j only. 

The individual charge-dynamical scalars are charge density, charge-weighted internal energy 
and pressure (per unit mass), 

P? = (59) 

a 

Eq = Y^^Ea = Yj^na-na)qac' = l^^aC^ / 7(7- 1) (60) 
a a a •'^u 

Pq = Y.—Pa = l^L^ai f{yy)^adClu (61) 

^ma 2>^ Ja^ 
and an enhancement jq in the spatial electric current due to relativistic beaming effects 

j' = Y^qa L y^y^adClu = = Jq^^ana^a = Iqj (62) 

a •'^u a a 

Because the partial internal energies and pressures pa have been defined previously (equations 
48 and 50), and because qa and are known, the quantities and Pq (equations 60 and 61) are not 
new variables but rather different weightings of known equations of state. Only the six quantities 
p^, j', and yq are new ones that need their own charge-dynamical equations. Five of those equations 
are, respectively, (56), the three components of (57) orthogonal to U, and the component of (57) 
parallel to U. Because j' is a four- vector, the sixth equation is a constraint on its components 

U-j' = (63) 
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similar to equation (12). 

The sources and sinks of current on the right-hand side of equation (57) are the Lorentz effect, 
the Hall effect, and resistive losses, with coefficients 



i = Y.i^lna/ma) (64) 

a 




where |j| = (— j • j)^^^ is the magnitude of the spatial current. 

Note that the definition of Jq (equation 62) makes use of the fact that \'a and are essentially 
parallel, resulting in the current j being enhanced by an average Lorentz factor jq (which is > 1). 
We can derive an equation for jq in a manner similar to that used for equation (55), arriving at 

j'-V.(Y,) = I^Tip, - l[j'-(U-V,U) +U-V,8,+ (e, + p,)V,-U] 

-Y,(Y,-l)Vx-j'- ^j'-F-U (67) 

While analogous to equation (55), equation (67) is of a very different character. Everything in it, 
including j', p^, U, F, and even Eq and pq can be considered known. What remains is a current- 
weighted gradient of the Lorentz factor equaling a quadratic function of that Lorentz factor. In the 
rest frame of the fluid (or when the fluid is at rest), the gradient j' • V^^ loses all time dependence, and 
the equation takes on an elliptical character. The distribution of Jq must be solved implicitly on 
each hypersurface, employing appropriate boundary conditions, etc. Equation (67) can be thought 
of as a constraint on the current beaming factor jq, enforcing the conservation of charge- weighted 
internal energy flow via j' at the same time as conservation of charge flow via j is enforced by 
equation (56). 

Implicit equations for Lorentz factors are not unusual in relativistic hydrodynamics or magneto- 
hydrodynamics (Duncan & Hughes 1994; Marti & Mueller 1999; Hughes et al. 2002; Koide 2003). 
In most cases, however, they are simple algebraic equations that need to be solved in a single cell 
in spacetime at each time step. In this case, the equation contains gradient information on jq and, 
therefore, must be solved over the entire hypersurface simultaneously. When the fluid is not at rest 
with respect to the frame in which the gradient Vjc is computed, and is flowing relativistically, j' 
can have a significant temporal component, rendering equation (67) an hyperbolic equation. Nev- 
ertheless, computationally, it would be wise to solve this particular equation implicitly at all time 
steps in order to avoid numerical problems when the local fluid velocity slows down. 
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5.3. Charge Dynamics in Component Form 



The equations of charge dynamics (56-58) are valid in any frame. Therefore, we can imme- 
diately convert them to component form in any metric. Using p^, j', and as the variables, they 
are 



[(p,U^+r/Y,)v^],, = 

j'^U^g^v = 
with the charge-current-pressure tensor given by 

C«P = [p,+ (e,+p,)/c2]u«uP + U«j'^+ruP + p,g«P 
and the usual Einstein summation convention and comma derivative applying 



(68) 
^ (69) 
(70) 

(71) 



/i=0 

In addition, we can replace the zeroth (temporal) component of equation (69) with that pro- 
jected along the four- velocity to get a component form of equation (67) for the current beaming 
factor 



1 



i^q+Pq) 



"U V-8),n -3-J ^/ivU 



-g 



(72) 



6. Discussion 



6.1. Special Cases 



It is useful to examine a few special cases of relativistic charge dynamics to compare with 
previous authors' work. 
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6.1.1. Steady State with No Hall Term 

Under many conditions the terms on the left hand side of (57) are small compared to those 
on the right hand side, because the time and length scales over which plasma properties vary are 
long. In addition, the Hall effect is often small compared to the Lorentz and resistive effects. 
The remaining terms, when projected orthogonal to U, then reduce to equation (14). They can 
be reduced further to the time-independent relativistic and nonrelativistic forms in section 2 under 
conditions of infinite conductivity, sub-relativistic flow, etc. 

6.1.2. Cold Plasma and Nonrelativistic Flow 

Equation (57) also can be reduced to the form (17) if we make the following cold plasma 
assumptions: 1) the charge-weighted internal energy and pressure are negligible compared with 
p^c^ and 2) the current beaming factor = 1 so that j' = j. These conditions allow the charge- 
current-pressure tensor to be re-written as 

C = + - Pq\}®\} + pqg 

and allow us to ignore the temporal component of (57), which is now redundant with the conser- 
vation of charge equation. With the assumption that the metric is that of Minkowskian flat space 
(with non-zero components g"" =[— 1) !» l],no sum on a), the three spatial components of 
(57) (those projected orthogonal to e?, not eu) to reduce to equation (17). 

If we additionally make the following nonrelativistic assumptions that U = [c, Vx, Vy, V^] and 
|y| << c, so that the orthogonal current j is approximately the spatial current / = [0, J^-, Jy, Jz], 
then equation (57) then reduces to equation (16). 



6.1.3. Relativistic Pair Plasma 

For highly relativistic flows a time-dependent form of equation (57) will be needed. However, 
for certain plasmas, such as a relativistic pair plasma near a black hole, equation (67) still can 
be simplified somewhat. In this case, m_ = m+ and q- = —q+, so Eg = pq = 0, and the Hall 
coefficient vanishes explicitly. Then equation (67) reduces to 

j' • V.(Y,) = 7^ [ P, - j' • (U • V,U)/c2 ] - Y,(Y, - 1) V, • j' (73) 

In other words, the beaming factor distribution is determined by the competition between local 
'creation' of charge by collisions and the loss of charge through the beamed current j'. 
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6.1.4. Uniform Adiabatic Index and Mixture 

More generally, if the plasma is made up of particle partial fluids that have the same adiabatic 
index, and if that index and the fractional pressure Ha = Pa/ P of each species are uniform through- 
out Q.X, then we have each quantity E^, £, and Eq related to their respective pressures by the simple 
relation 

ei = f^Pi (74) 

where i = a, q, or blank and T is the adiabatic index. Then the ratio of charge-weighted to total 
pressure (and that for internal energy) is a uniform constant throughout D.x 

P, ^ ^ = (75) 

P ^ a 

We then can use equation (55) to eliminate the charge-weighted thermodynamic quantities in (67) 
to obtain a simpler equation for 

j' • V.(Y,) = 7^ [ P, - (j' - CH/c^) . (U ■ V,U)/c2 - V, ■ (j' - CH/c2) ] 

+ Y,[V.-j'-^^^j'-F-U] (76) 

The pair plasma is a special case of these conditions, and equation (76) reduces to (73) in this case 
(^, h 0). Equation (76) is useful for showing how the equation for jq is decoupled from the 
terms involving Eq in equation (67). The quantity j' — ^H/c^ is a residual (beamed) current whose 
properties can affect the value of the current beaming factor jq. 



6.2. One-fluid vs. Multi-fluid Theory 

The advantage of one-fluid theory, of course, is that, by summing the multi-fluid equations 
and deriving thermodynamic quantities that close the system, the many equations for each particle 
species are reduced to only five (45 and 46). However, with the introduction of a set of five new 
one-fluid equations (56 and 57), one must ask whether it is still useful to deal with the one-fluid 
equations rather than the more instructive multi-fluid equations, particularly if there are only two 
fluids (ions and electrons). 

The answer still appears to be a qualified 'yes', although it is quite likely that the multi-fluid 
equations will become more important in the next few decades, if not sooner. First of all, most cur- 
rent astrophysical codes (hydrodynamic and magnetohydrodynamic, mainly) are one-fluid codes. 
While nevertheless a significant task, the addition of the one-fluid charge dynamical equations to 
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existing MHD codes is still much less effort than developing an entirely new multi-fluid code. 

Secondly, even when constructing a new code that involves only ions and electrons, developing 
that two-fluid code will be a much greater task than developing a one-fluid code with the ten equa- 
tions of hydro- and charge dynamics. In one-fluid theory, the collision terms can be treated with 
much less rigor than in a two-fluid code. In the former, because thermodynamic equilibrium is 
assumed, one need only postulate an approximate resistivity, as we have done in section 5. In the 
latter case, however, one must carefully handle coUisional momentum and energy transfer between 
each species, as well as the scattering of particles by electromagnetic oscillations with wavelengths 
shorter than a cell size. Otherwise, the simulated multi-temperature structure of the fluid will be 
meaningless. Finally, as is the case in simulations of late stages of stellar evolution and collapse, 
when dealing with the collapse of dense matter to black holes, there probably will be many more 
than two species of particle (neutrons, protons, electrons, positrons, heavy iron-peak nuclei, etc.). 
With five multi-fluid equations for each species, the number of equations to integrate could be 
significantly greater than the ten needed for hydro- and charge dynamics. 

Therefore, any multi-fluid astrophysical codes that are to be developed in the next few years 
are likely to be two-fluid only and probably will make the assumption of thermodynamic equi- 
librium initially. In that case their results will be similar to those obtained by older MHD codes 
that have been modified to handle charge dynamics. Nevertheless, these new codes will become 
increasingly sophisticated as more particle physics is added and should lead to a greater under- 
standing of black hole formation than is possible from one-fluid simulations alone. 

7. Conclusions 

This paper has used geometric frame-independent tensor notation to derive what the author 
believes is the first set of one-fluid equations to be published that correctly describes the flow of 
charge in general relativistic environments. Previously used or suggested approximations (e.g., 
force-free field, magnetohydrodynamics, even currently-available relativistic forms of the gener- 
alized Ohm's law) do not correctly describe black hole astrophysics on time scales much shorter 
than the collapse time (<< GM/c^) or in strong gravitational fields. The principle equations of 
charge dynamics are (56) and (57), with an alternative form for the relativistic current equation 
(67). Proper handling of charge flow in such environments will be important for understanding the 
details of highly relativistic astrophysical events like black hole formation and relativistic jet gen- 
eration, which may be important for understanding gamma-ray bursts, etc. (Multi-fluid equations 
are also derived, but the collision terms are not treated with sufficient rigor in this paper to make 
them useful for detailed simulations at this time.) These equations of 'charge dynamics', also given 
in component form in (68, 69, and 72), are valid for any flow velocity (causal), in any reference 
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frame (covariant), and in any spacetime metric (general relativistic). They therefore are suitable 
not only for flow in stationary metrics like Schwarzschild or Kerr, but also for general numerical 
relativity calculations that include a fluid with an embedded electromagnetic field. The equations 
of charge dynamics were shown to reduce to a variety of prior "generalized Ohm's laws" in the 
appropriate cold-plasma, flat-space or nonrelativistic limits. 

In addition to the general relativistic nature of the equations, the principal difference between 
this paper and previous work is that it does not make the assumption of a cold plasma. Not only 
must one deal with quantities such as charge-weighted internal energy and pressure, one also must 
solve for the current beaming factor that distinguishes the beamed current j' from the unbeamed 
j. The equation for y^, generated by the subtraction of the energy equations for positive and nega- 
tive ions, is primarily implicit and global in character, providing a constraint on the flow of charge. 
This result is in sharp contrast to the assertion (Khanna 1998) that a one-fluid theory is only possi- 
ble for a cold plasma. The charge-current tensor approach produces a fourth "energy" equation for 
the beamed current that can take into account fast current drift velocities and hot plasma as well as 
fast bulk speeds. 

While of course valid in more benign environments, the charge-dynamical form of Ohm's law 
probably will be needed only in very violent environments such as electromagnetic, rotating black 

hole formation. These events can have significant fluid and metric shear that may affect charge 
flow on time scales shorter than the current rise time. However, it is precisely these environments 
that are believed to obtain in "coUapsars", which form in the centers of massive stars, and in 
mergers of neutron stars in close binary systems. And it is these systems that are believed to lead 
directly to the formation of electromagnetic jets and their associated observable events. Charge 
dynamics, therefore, may play a significant role in understanding the engine that generates the 
highly relativistic jets seen in gamma-ray bursts and other sources. 

The author is grateful for discussions with M. Miller. He is also grateful for a JPL Institutional 
Research and Development grant, and for the continued hospitality of the TAPIR group at Caltech. 
The research described in this paper was carried out at the Jet Propulsion Laboratory, California 
Institute of Technology, under contract to the National Aeronautics and Space Administration. 



A. Velocity Integrals in Three-Space 

The velocity integrals in sections 4 and 5 were cast as ones over the four-dimensional volume 
Clu. This causes each integral, even those for quantities that appear in nonrelativistic dynamics (na. 
Pa, ^a), to contain a Lorentz factor (equation 30) in its kernel. In this Appendix we show that each 
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can be converted to its more familiar, three-space form by confining the integration to take place 
on the mass hyperboloid only. 

We begin by defining a three-dimensional density 

as the distribution of particles in three-momentum p, where is the particle rest mass and dTp = 
m^dTu. This distribution is related to by a delta function that enforces the mass hyperboloid 

= mlfahi^-y'c/l) ^ mlfa^^^^^ (A2) 

u"/c 

(See equation 20.) 

The four-volume integral of any kernel weighted by y, now can be replaced by a three- 
integral over fa on the mass shell, with no y weighting: 

= / f-'^'P [.^d{u^-yc)K{x,u^p)du^ 

= / K{x,yc,p)fad^p (A3) 

Equations (36), (37), (48), and (50) now take on their familiar forms 

na = I fad^P 



1 



p 



e« = / {y-\)ma(p- fad^p 

Pa = ^ p-yfad^P 

3 Jr„ 



(Recall that in the local Lorentz frame of the fluid v = (0, v\ v^, v^).) However, V'^ in equation 
(41), which is used to construct the beamed current j', has no nonrelativistic analog (other than j 
itself when |v| « c), and therefore must always involve the Lorentz factor 



V'« = - / yy fad' p 
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